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ABSTRACT 

Long time coverage and high radial velocity precision have allowed for the discovery of additional 
objects in known planetary systems. Many of the extrasolar planets detected have highly eccentric 
orbits, which raises the question of how likely those systems are to host additional planets. We inves- 
tigate six systems which contain a very eccentric (e > 0.6) planet: HD 3651, HD 37605, HD 45350, 
HD 80606, HD 89744, and 16 Cyg B. We present updated radial- velocity observations and orbital solu- 
tions, search for additional planets, and perform test particle simulations to find regions of dynamical 
stability. The dynamical simulations show that short-period planets could exist in the HD 45350 and 
16 Cyg B systems, and we use the observational data to set tight detection limits, which rule out 
additional planets down to a few Neptune masses in the HD 3651, HD 45350, and 16 Cyg B systems. 
Subject headings: extrasolar planets - planetary dynamics - stars: planetary systems 



1. INTRODUCTION 

One surprising result that has come out of the more 
than 200 extrasolar planet discoveries to date is the wide 
range of eccentricities observed. Unlike our own Solar 
system, many of the extrasolar planets which are not 
tidally locked to their host stars have moderate eccen- 
tricities (e > 0.2), and 15 planets have high eccentricities 
(e > 0.6). These observations have spawned several the- 
ories as to the origin of highly eccentric extrasolar plan- 
ets. One such method, planet-planet scattering, occurs 
when multiple jovian planets form several astronomical 
units (AU) from the host star and then interact, leav- 
ing one in an eccentr ic orbit and often ejecting the other 

asio fc Fordlll996f ). This method has been proposed 
to ex plain the a,rchite cture of the v And planetary sys- 
tem (|Ford et al.ll2005f ). which contains a hot Jupiter as 
well as two jovi an planets in moderately eccentric or- 
bits. iLin &: Idal ([1997) suggested a merger scenario m 
which inner protoplanets perturb each other and merge 
to form a single massive, eccentric planet with e^O.3 
and a - 0.5 - 1 AU. 

Interactions with stellar companions are another pos- 
sible way to boost a planet's eccentricity. Of the 15 
stars hosting a planet with e > 0.6, six are also known 
to possess stellar-mass co mpanions in wide bin ary or- 
bits: HD 3651 (Mugraucr st all l2006t Luhma n et all 
[20071) . HD 20782 (Dcsidcra fc Barbieri"20 07|), HP 80606 , 
HD 89744 (jWilson et al. 2001; Mug raue7eral] _ |2004f ) . 
16 Cyg B, and HD 222582 (jR,aghavan'"eraI1 \2(M) . If 
the inclination angle between the planetary orbit and a 
stellar companion is large, the Kozai mechanism (jKozail 
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fl96l) can induce large-amplitude oscillations in the ec- 
centricity of the planet (e.g. Malmberg et al. 2006). 
These oscillations can be damped by general relativis- 
tic effects and by interaction with other planets, and 
hence are most effective in systems with a single planet 
in an orbit a ^ 1 AU from the host star (jTakeda &: Rasiol 
|2005() . The Kozai mechanism has been suggested to ex- 



plain the high eccentric ity of 16 Cyg Bb (Holman et al 
1997; Mazch et al.TgW) and HD 80606b (Wu fc Murray 
2003). Ha uscr fc Ma rcv (.1999) found the inclination of 
16 Cyg B orbiting the system barycenter to lie between 
100 and 160 degrees, where 90 degrees is an edge-on ori- 
entation. However, it is the difference in inclination be- 
tween the orbital planes of the planetary and stellar com- 
panion that is critical in determining the importance of 
the Kozai mechanism, and the inclination of the planet's 
orbit is generally not known for non-transiting systems. 

Of the 192 known planetary systems, 23 (12%) 
are multi-planet systems. Recent discoveries of ad- 
ditional obj ects in systems known to host at least 
one planet (lUdrv et al] l2007t [ Wittcnniv eT et all l2007t 
iRivera et al MgOOa IVogt et al.ri2005 : .McArthur et all 
iU l2004f ) suggest that multiple-planet 
Of particular interest are sys- 
planet and a low-mass "hot 
75732), GJ 876, ^i Arae 



2004; San tos et al 
systems are common 
tems which host a jovian 
Neptune," e.g. 55 Cnc (=HD 
(=HD 160691), Gl 777A (=HD 190360). 



Motivated by 

the discoveries of hot Neptunes in known planetary sys- 
tems, we have undertaken an intensive survey of selected 
single-planet systems to search for additional low-mass 
companions. Three of the planetary systems discussed in 
this paper (HD 3651, HD 80606, HD 89744) are part of 
this campaign. The excellent radial-velocity precision of 
the High Resolution Spectrograph on the Hobby-Eberly 
Telescope (HET), combined with queue-scheduling, al- 
low us to time the observations in such a way as to min- 
imize phase gaps in the orbit of the known planet, and 
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also to act quickly on potential new planet candidates. 
T he use of the HET in this manner is discussed further 
in ICochran et al.l ()2004l ) with regard to the discovery of 
HD 37605b. 

In this work, we aim to combine observational lim- 
its on additional planets in known planetary systems 
with dynamical constraints obtained by N-body simu- 
lations. The observations address the question: What 
additional planets are (or are not) in these systems? The 
dynamical simulations can answer the question: Where 
are additional planets possible? Section 2 describes the 
observations and the test particle simulations for six 
highly eccentric planetary systems: HD 3651, HD 37605, 
HD 45350, HD 80606, HD 89744, and 16 Cyg B. We 
have chosen these systems based on two criteria: (1) 
Each hosts a planet with e > 0.6, and (2) Each has been 
observed by the planet search programs at McDonald 
Observatory. In §3, we present and discuss the results 
of the updated orbital fits, dynamical simulations, and 
detection limit computations. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Radial- Velocity Observations 

Five of the six stars considered in this work have 
been observed with the McDonald Observatory 9.2 m 
Hobby-Eberly Telescope (HE T) using i ts High Reso- 
lution Spectrograph (HRS) fl998h . A fuU de- 
sc ription of the HET p lanet search program is given 
in ICochran et al.l ()2004f ). For 16 Cyg B, observations 
from McDonald Observatory were obtained only with the 
2.7 m Harlan J. Smith (HJS) telescope; the long-term 
planet search program o n this telescope is described in 
I Wittenmyer et aIT(|2006l ). All available published data 
on these systems were combined with our data from Mc- 
Donald Observatory in the orbit fitting procedures. 

2.2. Numerical Methods 

To place constraints on the architecture of planetary 
systems, we would like to know where additional objects 
can remain in stable orbits in the presence of the known 
planet (s). We performed test partic le simulations us- 
ing SWIFT^ ()Levison fc Duncanlll994l ) to investigate the 
dynamical possibility of additional low-mass planets in 
each of the six systems considered here. Low-mass plan- 
ets can be treated as test particles since the exchange 
of angular momentum with jovian planets is small. We 
chose the regularized mixed-variable symplectic integra- 
tor (RMVS3) version of SWIFT for its ability to handle 
close approaches between massless, non-interacting test 
particles and planets. Particles are removed if they are 
(1) closer than 1 Hill radius to the planet, (2) closer 
than 0.05 AU to the star, or (3) farther than 10 AU 
from the star. Since the purpose of these simulations 
is to determine the regions in which additional planets 
could remain in stable orbits, we set this outer bound- 
ary because the current repository of radial- velocity data 
cannot detect objects at such distances. 

The test particle simulations were set up following the 
methods used in Barnes & Raymond (2004), with the 
exception that only initially circular orbits are consid- 
ered in this work. For each planetary system, test par- 

^ SWIFT is publicly available at 

|http: //www. boulder. swri.edu/^hal/swift.htmlj 



tides were placed in initially circular orbits spaced ev- 
ery 0.002 AU in the region between 0.05-2.0 AU. We 
have chosen to focus on this region because the dura- 
tion of our high-precision HET data is currently only 
2-4 years for the objects in this study. The test parti- 
cles were coplanar with the existing planet, which had 
the effect of confining the simulation to two dimen- 
sions. Input physical parameters for the known planet 
in each system were obtained from our Keplerian or- 
bit fits describ ed in § 3.1, and from recent literature 
for 16 Cyg B (iWittenmver et al.l l200l and HD 45350 
(|Endl et al.ll2006( l. The planetary masses were taken 
to be their minimum values (sin i = 1). The systems 
wereintegrated for 10^ yr, following iBarnes fc Ravmondl 
(2004) and allowing completion of the computations in 
a reasonable time. We observed that nearly all of the 
test-particle removals occurred within the first 10^ yr; 
after this time, the simulations had essentially stabilized 
to their final configurations. 

3. RESULTS AND DISCUSSION 

3.1. Updated Keplerian Solutions for 4 Systems 

We present updated Keplerian orbital solutions for 
HD 3651b, HD 37605b, HD 80606b, and HD 89744b 
in Table 1. A summary of the data used in our anal- 
ysis is given in Table 2, and the HET radial veloci- 
ties are given in Tables 3-6. The velocity uncertain- 
ties given for the HET data represent internal errors 
only, and do not include any external sources of error 
such as stellar "jitter." The parameters for the remain- 
ing t wo planets, HD 4535 0b and 16 Cyg Bb, are taken 
from lEndl et~al] (|2OO60 and IWittenmver et all (|2007f ). re- 
spectively. Radial velocity measuremen ts from the HET 
are given for HD 45350 in lEndl et all (12006), and ve- 
locitie s for 16 Cyg B from the HJS telescop e are given 
in Wittenmver et al.l ()2007f ). As in Wittcn mver et ahl 
(|2007h . all available published data were combined with 
those from McDonald, and the known planet in each 
system was fit with a Keplerian orbit using GaussFit 
(jJeffervs et al.|[T987l ) , allowing the velocity offset between 
each data set to be a free parameter. Examination of the 
residuals to our Keplerian orbit fits revealed no evidence 
for additional objects in any of the six systems in this 
study. 

The Saturn-mass (M sini = 0.2M,T„p ) planet HD 3651b 
was discovered by Fis cher et al.l ()2003l ) using observa- 
tions from Lick and Keck. We fit these data, which were 
updated in lButler et al.l (j2006 D. in combination with ob- 
servations from the HJS and HET at McDonald Obser- 
vatory. The HET data, which consist of multiple expo- 
sures per visit, were binned using the inverse- variance 
weighted mean value of the velocities in each visit. The 
standard error of the mean was added in quadrature to 
the weighted rms about the mean velocity to generate 
the error bar of each binned point (N=29). The rms 
about the combined fit for each dataset is: Lick & Keck- 
6.6 m s-\ HET-9.4 m s^^, HJS-12.2 m s^^. The fitted 
orbital parameters for H P 3651b are o f comp arable pre- 
cision to those reported in lButler et al.l ()2006f ). and agree 
within 2cr. T he recent discovery of a T dwarf companion 
to HD 3651 (|Mugrauer et al.l 120061: ILuhman et al.ll2007D 
prompts an interesting exercise: Can the radial-velocity 
trend due to this object be detected in the residuals after 
removing the planet? We detect a slope of —0.27 ± 0.05 
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m s^^ yr^^, indicating that we are indeed able to discern 
a trend which is possibly due to the binary companion. 
However, the reduced of the orbital solution is not 
significantly improved by the inclusion of a linear trend 
(Ax^=0.f8). The parameters given in Table 1 were ob- 
tained from the fit which did not include a trend. 

We present 23 new HET observ ations for HP 37605 ob - 
tained since its announcement bv lCochran et al.l (|2004D . 
The data now span a total of 1065 days. The best fit 
is obtained by including an acceleration of —20.5 ± 2.1 
m s~^ yr~^, indicating a distant orbiting body. Such a 
finding would lend support to the hypothesis that very 
eccentric single-planet systems originate by interactions 
within a wide binary system. The shortest period that 
this outer companion could have and still remain consis- 
tent with the observed acceleration and its uncertainty 
over the timespan of the observations is about 40 yr, as- 
suming a circular orbit. This object would then have a 
minimum mass in the brown dwarf range. 

The planet or biting HD 80606, first announced by 
iNaef et al] ()2001h . is the most eccentric extrasolar planet 
known, with e = 0.933 ± 0.001 (Table 1). We have fit 
the CO RALIE data in comb ination with the Keck data 
given in iButler et al.l ()2006f ) and 23 observations from 
HET. The extreme velocity variations induced by this 
planet greatly increase the sensitivity of orbit fits to the 
weighting of individual measurements. Since the uncer- 
tainties of the HET velocities given in Tables 3-6 repre- 
sent internal errors only, we experimented with adding 
1-7 m s~"'^of radial-velocity "jitter" in quadrature before 
fitting the data for HD 80606. For all of these jitter 
values, the fitted parameters remained the same within 
their uncertainties. Table 1 gives the param eters derived 
from a fit which added 3.5 m s~^of jitter (jButler et al.l 
I2006D to the HET data. The rms about the combined 
fit is: CORALIE 18.7 m s^^, HET-7.5 m s^^, Keck-5.6 
m s~^. iButler et al.l (|2006l ) noted that the eccentricity e 
and the argument of periastron u! had to be held fixed 
in their fit to the Keck data alone. However, the large 
number of measurements included in this work allowed 
GaussFit to converge with all parameters free. 

For HD 89744b, we combine data from the HET with 
6 me asurements from th e HJS telescope and Lick data 
from IButler et al.l ()2006f ). The HET data were binned 
in the same manner as for HD 3651, resulting in N=33 
independent visits. The rms about the combined fit for 
each dataset is: Lick-17.1m s"\ HET-10.7m s"\ HJS~ 
9.5 m s~^. As with HD 3651b, o ur de rived parameters 
agree with those of IButler et al.l (|2006f) within 2a. The 
scatter about our fit remains large, most likely due to the 
star's early spectral type (F7V), which hinders precision 
radial- velocity measurements due to the smaller number 
of spectral lines. For example, the F7 y star HD 22128 7 
was recently found to host a planet (jNaef et al.ll2007f ): 
despite the superb instrumental precision of the HARPS 
spectrograph, that orbital solution has a residual rms of 
8.5 m s-i. 

3.2. Test Particle Simulations 

The results of the dynamical simulations are shown in 
Figures 1-3. The survival time of the test particles is plot- 
ted against their initial semimajor axis. As shown in Fig- 
ure 1, the short-period planets HD 3651 and HD 37605 
sweep clean the region inside of about 0.5 AU. In both 



of these systems, however, a small number of test parti- 
cles remained in low-eccentricity orbits near the known 
planet's apastron distance, near the 1:2 mean-motion res- 
onance (MMR). In the HD 3651 system, particles re- 
mained stable beyond about 0.6 AU, which is not sur- 
prising given the low mass of the planet. For HD 37605, 
two distinct strips of stability are seen in Fig. 1, corre- 
sponding to the 1:2 and 1:3 MMRs. The eccentricity of 
the test particles in the region of the 1:2 MMR oscillated 
between 0.00 and 0.06. Particles in 1:3 MMR oscillated 
in eccentricity with a larger range, up to e ~ 0.4, which 
is expected due to secular forcing. As with HD 3651, the 
region beyond about 0.8 AU was essentially unaffected 
by the planet. 

Figure 2 shows the results for the HD 45350 and 
HD 80606 systems. The long period (963.6 days) and 
relatively large mass (M sin i—1.8 Mjup) of HD 45350b 
restricted stable orbits to the innermost 0.2 AU. These 
test particles oscillated in eccentricity up to e ~ 0.25. 
The 4Mjup planet orbiting HD 80606 removed all test 
particles to a distance of about 1.5 AU, and only be- 
yond 1.75 AU did test particles remain in stable orbits 
for the duration of the simulation (10^ yr). A region of 
instability is evident at 1.9 AU due to the 8:1 MMR. 
Figure 3 shows that HD 89744b eliminated all test par- 
ticles except for a narrow region near the 8:3 resonance. 
For the 16 Cyg B system, only particles inside of about 
0.3 AU remained stable, leaving open the possibility of 
short-period planets. The surviving particles oscillated 
in eccentricity up to e 0.45, but these simulations treat 
the star as a point mass, and hence tidal damping of the 
eccentricity is not included. Our results are consistent 
with those of iMenou fc TabachnikI (|2003l ). who investi- 
gated dynamical stability in extrasolar planetary systems 
and found that no test particles survived in the habitable 
zones of the HD 80606, HD 89744, and 16 Cyg B systems. 

3.3. Detection Limits 

Three of these systems (HD 3651, HD 80606, 
HD 89744) were monitored intensely with the HET as 
part of a larger effort to search for low-mass, short pe- 
riod planets. No evidence was found for any such objects 
in these or any of the six systems in this work. We then 
asked what limits can be set on additional planets us- 
ing the high-precision HET data we have obtained. The 
procedure for determining cornpanion limits was identi- 
cal to the method described in iWittenmver et^Il (|2006f) . 
except that in this work, the best-fit Keplerian orbit for 
the known planet (see § 3.1) was removed before perform- 
ing the limits computations. In this way, we determined 
the radial- velocity amplitude K for which 99% of planets 
would have been detected in the residuals. The eccentric- 
ity of the injected test signals was chosen to be the mean 
eccentricity of the surviving particles from the simula- 
tions described in § 3.2. Only the regions in which test 
particles survived were considered in these limits compu- 
tations. 

The results of these computations were highly varied, 
reflecting the differing observing strategies employed for 
these six objects. In particular, HD 3651, HD 80606, 
and HD 89744 were monitored intensely with the HET 
as part of a search for short-period objects, whereas 
HD 37605 and HD 45350 were only observed sporadically 
after the known planet orbits were defined and published 
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(iCochran et alJ [200l lEndl et alJ 120061) . and 16 Cyg B 
has only been observed with the HJS telescope at a fre- 
quency of at most once per month. The companion limits 
are shown in Figures 4-6; planets with masses above the 
solid line can be ruled out by the data with 99% confi- 
dence. Not surprisingly, the tightest limits were obtained 
for HD 3651 (Figure 4), which had a total of 195 mea- 
surements, including 29 independent HET visits. For 
periods less than about 1 year, we can exclude planets 
with M sin i ^ 2 Neptune masses. Similar results were 
obtained for 16 Cyg B (N=161), where the limits ap- 
proach a Neptune mass (Figure 6). Since the detection 
limits generally improve with the addition of more data 
and with higher-quality data, we can define a quantity 
to measure the goodness of the limits. A simple choice 
would be iV/tJ, where N is the total number of observa- 
tions, and a is the mean uncertainty of the radial- velocity 
measurements. The values of N and tr are given in Ta- 
ble 2. 

In the HD 45350 system, the results of the dynamical 
simulations complement those of the detection limit de- 
terminations. Very tight limits are obtained in close or- 
bits (a ^ 0.2 AU). In this region, test particles were stable 
(Fig. 2) and our observations can exclude planets with 
M sin i between about 1 and 4 Neptune masses. Similar 
results were obtained for the 16 Cyg B system, in which 
test particles remained stable inward of a ~ 0.3 AU. In 
that region, planets of 1-3 Neptune masses can be ex- 
cluded by our limits determinations (Fig. 6). In most 
of the limits determinations, there are multiple "blind 
spots" evident where the periodogram method failed to 
significantly recover the injected signals. Typically this 
occurs at certain trial periods for which the phase cov- 
erage of the observational data is poor, and often at the 
1-month and 1-year windows. 

For none of HD 37605 (Fig. 4), HD 80606 (Fig. 5), or 
HD 89744 (Fig. 6) could additional companions be ruled 
out below about 0.7 Mjup, and for most orbital periods 
tested, the limits were substantially worse. One possible 
explanation for this result is that the sampling of the 
observations was poorly distributed in phase for many 
of the injected test signals, making significant recovery 
by the periodogram method difficult. This is evidenced 
by the "jagged" regions in the plots. Also, the intrinsic 
scatter for those three systems was too large to permit 
tight limits determination. This is certainly reasonable 
for the F7 star HD 89744. The three systems with the 
best hmits (HD 3651, HD 45350, and 16 Cyg B) also 
had the lowest rms scatter about their orbital solutions 
(mean=8.9 ± 1.4 m s~^; Table 1). In contrast, the mean 
rms for the remaining three systems was 13.7±0.6 m s^^. 
Additional factors such as a paucity of data (HD 37605) 
and short time baselines (HD 80606, HD 89744) made the 
determination of useful companion limits challenging for 
some of the planetary systems in this study. 

4. SUMMARY 

We have shown that for a sample of six highly eccen- 
tric extrasolar planetary systems, there is no evidence for 
additional planets. Test particle simulations show that 
there are regions detectable by current surveys (i.e. for 
a < 2 AU) where additional objects can exist. For 
HD 3651 and HD 37605, we find that protected reso- 
nances are also present. Combining these simulations 



with detection limits computed using new high-precision 
HET data combined with all available published data is 
particularly effective for the HD 3651 and HD 45350 sys- 
tems. Additional short-period planets can be ruled out 
down to a few Neptune masses in the dynamically stable 
regions in these systems. 
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TABLE 1 
Keplerian Orbital Solutions 



Planet Period 


To 


e 




K 


M sin i 


a rms 


(days) 


(JD-2400000) 




(degrees) 


(in s^^) 


(Mj,p) 


(AU) m s-i 


HD 3651 b 62.197±0.012 


53932.2±0.4 


0.630±0.046 


250.7±6.3 


15.6±1.1 


0.20±0.01 


0.280±0.006 7.1 


HD 37605 b 55.027±0.009 


52992.8±0.1 


0.677±0.009 


218.4±1.7 


201.5±3.9 


2.39±0.12 


0.263±0.006 13.0 


HD 45350 b 963.6±3.4 


51825.3±7.1 


0.778±0.009 


343.4±2.3 


58.0±1.7 


1.79±0.14 


1.92±0.07 9.1 


HD 80606 b 111.428±0.002 


53421. 928±0.004 


0.933±0.001 


300.4±0.3 


470.2±2.5 


4.10±0.12 


0.460±0.007 13.5 


HD 89744 b 256.78±0.05 


53816.1±0.3 


0.689±0.006 


194.1±0.6 


263.2±3.9 


7.92±0.23 


0.91±0.01 14.4 


16 Cvp- R h 799 5+0 6 


50539.3±1.6 


0.689±0.011 


83.4±2.1 


51.2±1.1 


1.68±0.07 


1 fis+n 03 10 6 






TABLE 2 










Summary of Radial- Velocity Data 






Star 


N 


(T (m s 




AT (days) 


1 rT*t^ 


HD 3651 


163 


3.4 








Butler et al. f2006l 


HD 3651 


3 


6.1 








HJS^ 


HD 3651 


29 


2.1 








HET'' 


HD 3651 (total) 


195 


3.2 




7083 






HD 37605 (total) 


43 


2.9 




1065 




HET 


HD 45350 


38 


2.8 








Butler et al. f2006) 


HD 45350 


28 


4.2 








HET 


HD 45350 


47 


8.9 








HJS 


HD 45350 (total) 


113 


5.7 




2265 






HD 80606 


61 


13.7 








Naef et al. (2001) 


HD 80606 


46 


5.1 








Butler et aLlWOG) 


HD 80606 


23 


2.5 








HET 


HD 80606 (total) 


130 


8.7 




2893 






HD 89744 


50 


11.2 








Butler et al. (2006) 


HD 89744 


33 


3.2 








HET 


HD 89744 


6 


9.4 








HJS 


HD 89744 (total) 


89 


8.1 




2687 






16 Cyg B 


95 


6.3 








Butler et al. (2006) 


16 Cyg B 


29 


19.7 








HJS Phase 11'= 


16 Cyg B 


37 


7.4 








HJS Phase HI 


16 Cyg B (total) 


161 


9.0 




6950 







^ McDonald Observatory 2.7 m Harlan J. Smith Telescope.^ McDonald Observatory 9.2 m Hobby-Ebcrly Telescope.*^ Phase II indicates an earlier 
instrument setup detailed in lWittenmver et al Phase III is the current configuration. 
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Fk;. 1. Left panel: Survival time as a function of initial semimajor axis for test particles in the HD 3651 system after 10^ yr. The 
filled regions indicate test particles which survived. The orbital excursion of HD 3561b is indicated by the horizontal error bars at the top. 
Particles were placed on initially circular orbits with 0.05 < a < 2.00 AU. For all systems, the known planet removed particles which crossed 
its orbit. The dark region near 0.5 AU shows the stable 1:2 mean-motion resonance (MMR). Right panel: Same, but for the HD 37605 
system. The dark regions near 0.45 AU and 0.6 AU show the stable 1:2 and 1:3 MMRs. 
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Fig. 2.— Same as Fig. 1, but for the HD 45350 (left) and HD 80606 (right) systems. 
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Fig. 3.— Same as Fig. 1, but for the HD 89744 (left) and 16 Cyg B (right) systems. 
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Fig. 4. — Left panel: Detection limits for additional planets in orbits with e = 0.18 in the HD 3651 system. Planets in the parameter 
space above the plotted points are excluded at the 99% confidence level. The horizontal dashed line indicates the mass of Neptune. Right 
panel: Same, but for planets with e = 0.20 in the HD 37605 system. 
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Fig. 5. — Left panel: Detection limits for additional planets with e = 0.22 in the HD 45350 system. The horizontal dashed line indicates 
the mass of Neptune. Right panel: Same, but for planets with e = 0.31 in the HD 80606 system. Planets in the parameter spa<;e above the 
plotted points are excluded at the 99% confidence level. 
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Fig. 6. — Left panel: Detection limits for additional planets with e = 0.28 in the HD 89744 system. Right panel: Same, but for planets 
with e = 0.23 in the 16 Cyg B system. The horizontal dashed line indicates the mass of Neptune. Planets in the parameter space above 
the plotted points are excluded at the 99% confidence level. 
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TABLE 3 

HET Radial Velocities for HD 3651 



JD-2400000 Velocity (m s~^) Uncertainty (m s~l) 



53581.87326 


-19.1 


2.9 


53581 87586 


-19.4 


2.7 


5S581 87846 


-20.7 


2.7 


5Sfinn 7Q6fiQ 


-11.5 


2.4 


53600 79860 


-15.5 


3.0 


5S600 80050 


-22.8 


2.9 


5S604 79166 


-15.8 


1.9 


5S604 79356 


-18.8 


2.1 


5S604 79548 


-21.3 


2.1 


5S606 78169 


-19.3 


1.8 


5S606 7SS60 


-14.8 


2.1 


5S606 78551 


-24.0 


1.8 


53608 77236 


-18.8 


1.9 


53608 77426 


-18.0 


1.9 


53608 77617 


-18.8 


1.8 


53615 96280 


-28.0 


2.6 


53615 96471 


-31.9 


2.4 


5361 5 96662 


-37.8 


2.5 


53628 74050 


-6.8 


2.2 


53628 74240 


-14.5 


2.4 


53628 74431 


-5.5 


2.2 


53669 61012 


-18.2 


2.1 


53669 61 203 


-19.2 


2.2 


53669.61394 


-17.7 


2.4 


53678 78954 


-10.6 


2.4 


53678 79141 

rJ^JU c \J • 1 Cirri 


-8.6 


2.3 


53678 79332 


-2.3 


2.1 


53682 78423 


-15.4 


2.2 


53682.78609 


-15.0 


2.3 


53682 78801 


-11.9 


2.3 


53687 77684 


11.3 


2.2 


53687 77875 


8.7 


2.2 


53687 78066 


15.9 


2.2 


53691.75967 


9.6 


2.2 


53691 76158 


20.3 


2.1 


53691 76349 


15.9 


2.0 


53696 75837 


16.1 


1.8 


53696 76028 


18.6 


1.8 


53696 76220 


20.0 


2.0 


53694 75275 


18.0 


1.9 


53694 75466 


15.1 


2.0 


53694 75656 


17.8 


2.0 


53955 83401 


-0.5 


1.9 


53955 83593 


-1.2 


2.0 


53955 83785 


1.3 


1.9 


53956 82850 


0.4 


2.0 


53956 83046 


-1.0 


2.0 


,53956 832.36 


-5.4 


2.2 


53957 82201 


-2.1 


2.0 


53957 82392 


-1.3 


2.0 


53957 82583 


-3.6 


2.0 


53973 80721 


9.8 


7.3 


53973.81020 


3.5 


2.3 


,5.3973 81200 


-3.5 


2.0 


53976 78393 


-10.4 


2.4 


53976 78586 


-5.4 


2.1 


53976 78778 


-6.7 


2.3 


53978.97197 


-3.8 


2.6 


53985 95886 


-9.0 


2.3 


53985 96079 


4.3 


3.3 


53987 95335 


-8.3 


2.2 


53987 95527 


-8.0 


2.2 


53987.95719 


-12.0 


2.3 


53989 7381 7 


-13.2 


2.2 


53989.74009 


-13.2 


2.1 


53989.74203 


-18.6 


2.1 


54003.70719 


2.0 


2.2 


54003.70915 


4.7 


2.4 


54005.68297 


7.0 


2.5 


54005.68488 


11.1 


2.0 


54005.68690 


10.2 


2.1 


54056.77919 


-7.5 


2.2 


54056.78110 


-11.5 


2.1 


54056.78302 


-9.6 


2.3 


54062.55119 


20.1 


1.8 


54062.55312 


21.9 


2.0 


54062.55505 


20.9 


2.0 


54064.54710 


12.8 


2.0 


54064., 54902 


16.7 


2.1 
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TABLE 4 

HET Radial Velocities for HD 37605 



JD-2400000 


Velocity (m s ) 


Uncertainty (m s ) 


ci'lnno fi7l "SI 


487.6 


3.8 




495.5 


3.0 


(JOUUU . UUZrUO 


496.2 


3.0 


tJOUUO.UU4:U ( 


OUl .o 


9 Q 




499.8 


2.9 




482.1 


2.6 


c;'?n49 727Q7 


269.7 


2.8 




489.0 


2.6 


OOUUO . IJ'4iUo4: 


A7Q n 




OOU 1 L.U'-tOOO 


463.8 


2.6 




460.4 


2.6 




422.8 


2.5 


"jouoo-oyoou 


422.2 


2.8 


'^'^OSK '^Q'?7J^ 
oouoo.oyo 1 o 


418.6 


4.0 


■^mnSQ '=lQ'=17fi 

tJOUOi/.Ut/O 1 u 


379.1 


2.2 




343.7 


2.5 


uouyi.uoutjo 


323.2 


2.4 




302.1 


2.4 




302.1 


3.2 




193.8 


2.7 




164.9 


3.0 


'^'i9fif^ QA7AA 


112 9 


o.u 




113.2 


3.7 




74.6 


3.6 


'^'^9fifi Q7'^Qfi 
oo^uu.y i oyu 


119.2 


8.0 


'^'^9^^'^ Q99zl1 


471.6 


2.7 


'^'^'^ 1 S ^^1 Q97 


213.3 


3.0 


f^QQQp; Q91 SI 


496.9 


2.6 




493.9 


2.6 


'^'^'^77 ^^1 Qzll 


109.1 


2.7 


'^'^'^7S ^^1 1 SQ 
ooo ( o.o±xoy 


214.6 


2.7 




338.3 


2.6 


53381.64429 


436.1 


2.7 


53384.64654 


482.9 


2.8 


53724.85584 


468.2 


2.6 


53731.69723 


435.4 


2.7 


53738.67472 


404.3 


2.6 


53743.81020 


400.5 


2.6 


53748.64724 


348.4 


2.7 


54039.85015 


272.5 


3.1 


54054.96457 


437.4 


2.7 


54055.95279 


422.0 


2.9 


54067.76282 


376.4 


2.6 
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TABLE 5 

HET Radial Velocities for HD 80606 



JD-2400000 


Velocity (m s ) 


Uncertainty (m s ) 




-20.8 


3.0 


ooooo . u^uoy 


-49.5 


2.7 


^'H'^'SQ 82400 


-60.4 


3.0 


oooui.uzyoo 


f\A 7 




(joouo.uou 1 y 


-77 A 


2.4 


ooo / o.yo^ozr 


-88.4 


3.0 


ci'^Qyy sfi! 1 9 

'J'-JiJ 1 1 .OUiJ-i; 


-105.5 


2.4 


ooo 1 y. i ozrou 


-109.3 


2.7 


'^'^'^80 7A1 70 
ooooy . ( 'ii 1 u 


--L lU.O 




'^'^'^Ql 74400 
oooy ± . / 'iiuu 


-129.4 


2.4 


c;QQqf^ 797R'? 


-146.4 


2.3 


c;qqqq 791^1 ft 
oooyy. 1 ^o±o 


-158.4 


2.5 


53401.72497 


-174.7 


2.7 


53414.67819 


-219.8 


3.0 


53421.85529 


261.0 


2.2 


53423.86650 


322.1 


2.0 


53424.85231 


245.9 


2.1 


53432.87120 


87.5 


1.9 


53433.60628 


70.0 


2.1 


53446.79322 


4.5 


1.9 


54161.85400 


-109.5 


2.8 


54166.83797 


-119.3 


2.4 


54186.76189 


-184.2 


2.3 
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TABLE 6 

HET Radial Velocities for HD 89744 



JD-2400000 


Velocity (m s ) 


Uncertainty (m s ) 




-184.5 


2.3 


(JO 1 ZO.OOJ-OO 


-238.6 


2.2 


(JO ( ZO.O(JOU 1 


-238.2 


2.5 


(JO I ZO.OUO'iU 


997 7 


9 

z.o 


(JO 1 Zr J .o'-tOi/^ 


-238.9 


2.5 


(JO I Z I .O'itJ I o 


-244.9 


2.4 


(JO 1 Z 1 . O'l / (JZr 


-242.9 


2.6 


(JO ( OU.OiOO 1 


-257.6 


2.5 


(JO / ou.oz±uu 


-248 2 


9 Q 
z.y 


(JO I OU.OZO±!j 


-253.4 


2.4 


<JO 1 OO.VJOiiUi 


-246.7 


2.8 


OO 1 OO.UOI^X 


-243.3 


2.4 


OO 1 oo.uouzu 


-236.0 


2.5 


OO 1 oo.ouoou 


940 
-z^U.o 


9 R 
z.u 


OO I OO.O-LU'iU 


-258.9 


2.4 


OO 1 00.0-Lt;J.i7 


-249.3 


2.5 


OO I 01.0-L 1 £70 


-242.8 


2.6 


'=m7^4 81 Q7^ 
OO 1 o'i.oj.y 1 o 


-243.9 


2.8 


'iS7S4 821 ^2 


-248.5 


2.4 


'^'^749 7Q1 1 Q 


-252.0 


2.8 


^'\7A9 7Q9QQ 
OO i iz. ( yzryy 


-257.2 


2.8 


(JO / ( y4: / y 


9'?Q 7 
-zoy. ( 


9 8 
z.o 


'^'i7'M 781 QQ 
ooj o i-. 1 oxyy 


-257.4 


2.9 


'^'i7f^1 7S'i78 

OO 1 o± . 1 OO 1 o 


-263.1 


2.5 


OO i o±. 1 oooo 


-268.0 


2.3 


OO i OO. t Oi-OO 


-273.1 


2.5 


^'\7^'^ 78^^81 

OO i OO. 1 OOOl. 


-278.7 


2.5 


'^'^7'^'^ 78(^07 

OO 1 OO . 1 0\)U 1 


-266.4 


2.4 


OO 1 0(J. 1 UVJOO 


-286.6 


2.3 


'^'^7'^'^ 7fi91 8 
(JO 1 OO . 1 uz±o 


-266.5 


2.6 


^'^7^^ 7fi'^Q7 

OO 1 OO . ( OOt? ( 


-274.9 


2.7 


OO 1 'xU.OJ-OULI 


-257.1 


1.9 


■=1^74^ 81 778 

OO I ^U.O-L 1 1 O 


-250.9 


2.1 


OO I ftyJ-O^UUX 


-245.2 


2.3 


f^^7^7 77nn9 

OO I O 1 • 1 I uu^ 


-277.6 


2.4 


OO 1 O 1 • 1 f xox 


-280.3 


2.4 


OO t O t • i 1 ouu 


-288.7 


2.2 


OO ( y I . u'iuuy 


-439.8 


3.1 


OO 1 cf 1 .vJrtOOT: 


-462.6 


2.8 


OO 1 y 1 .uouoy 


-452.5 


2.9 


'^'^SnQ fi9498 


-uoo.u 


2 4 


(joouy.uz i uu 


-uoo.o 


9 ^ 
z.o 


'^'i8nQ fi9Q79 
ooouy.uzy i z 


-659.2 


2.3 




-304.3 


3.0 


!jiJO'J 1 . 1 UO 1 w 


-324.0 


2.9 


'^'^8'^7 787*^1 
OOOO 1 . 1 O 1 1. 


-308.6 


2.7 


^'^^'^7 70077 
OOOO 1 . ( yu 1 ( 


-285 2 


9 R 

z.u 


'^'i8fifi fiQQ87 
ooouu.uyyo / 


-215.9 


1.7 




-228.3 


1.7 


OOOUyJ. 1 UU 1 U 


-220.4 


1.8 


ooouo.uoo'iy 


-251.6 


3.8 


ooouo.uoou^ 


-208.6 


2.9 


■iSSfiS fi8777 

OOOUO.UO 1 1 1 


-247.4 


9.7 


ooo ( o.uuyou 


-215.7 


1.6 


(JOOOO . uoouo 


-213.8 


1.8 


(JOOOO. uooo t 


-209.2 


1.7 


ooooo.uu-Luy 


-200.4 


1.7 


oooyu.uo / 1 u 


-203.4 


1.7 


oooyu . uoyoi 


-202.6 


1.9 


■^'^800 R4 1 '^4 
oooyu.U'i-Loi 


-203.2 


1.9 


oooyo.uzyoy 


-193.8 


2.0 


OOOifO.\)0\-0<y 


-189.3 


1.9 


OOOC/iJ.OOOJ-O 


-189.7 


1.8 


•^4047 Q481 1 

OrlUrl: 1 . yrtO J. i- 


-375.2 


4.8 


^4047 Q4QQ1 


-353.2 


4.5 


•^4047 Q'^1 79 
otu^ 1 . y o X ( z 


-362.6 


4.4 


'^40 "^0 QR948 
o^UoU . y oz^o 


41 n 


9 R 
z.u 


o^uou.yu^uo 


-423.0 


2.5 


54050.96657 


-420.1 


2.4 


54052.96488 


-426.8 


2.3 


54052.96762 


-437.1 


2.5 


54052.97035 


-447.6 


2.5 


54056.94606 


-468.0 


3.0 


54056.94786 


-466.4 


2.6 


54056.94964 


-479.4 


2.8 


54063.92981 


-599.1 


2.1 


54063.93166 


-594.8 


2.3 


54063.93348 


-592.3 


2.4 


54073.91213 


-685.8 


2.8 



